Global versus local billiard level dynamics: The limits of universality 
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Level dynamics measurements have been performed in a Sinai microwave billiard as a function 
of a single length, as well as in rectangular billiards with randomly distributed disks as a function 
of the position of one disk. In the first case the field distribution is changed globally, and velocity 
distributions and autocorrelation functions are well described by universal functions derived by 
Simons and Altshuler. In the second case the field distribution is changed locally. Here another 
type of universal correlations is observed. It can be derived under the assumption that chaotic wave 
functions may be described by a random superposition of plane waves. 
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In 1972 Edwards and Thouless noticed that the con- 
ductivity of a disordered system is closely related to the 
sensitivity of its eigenvalues on an external perturbation 



For a ring with a perpendicularly applied magnetic 
field they conjectured that the conductivity C is pro- 
portional to the averaged curvature of the eigenvalues, 

^ ~ (l%^H<p=o^! where tp is the magnetic flux through 
the ring. In 1992 Akkermans and Montambaux showed 
that the conductivity may alternatively be expressed in 

terms of the eigenvalue velocities, C ~ (^^§^\ 2 ^j @- This 
suggests to rescale the parameter and the eigenvalues by 
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where A is the mean level spacing. Szafer, Simons and 
Altshuler studied a number of parametric correlations of 
the rescaled eigenenergies [|]|| , in particular the velocity 
autocorrelation function 
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originally introduced by Yang and Burgdorfer [||, and 
conjectured a universal behavior as long as the so-called 
zero-mode approximation holds, i.e., in the range where 
the energy fluctuations show random matrix behavior. 
For the velocity distribution Simons and Altshuler found 
a Gaussian behavior || . The same behavior has been ob- 
tained by a completely different approach starting from 
the analogy between the level dynamics of a chaotic sys- 
tem and the dynamics of a one-dimensional gas with re- 
pulsive interaction . In the region of onset of local- 
ization deviations from the Gaussian behavior are found 

i-. 

Since in the zero-mode approximation the energy cor- 
relations of a disordered system are identical to that of 
random matrices, it came as no surprise that the uni- 
versal behavior of parametric correlations was found in 
billiard systems as well ||. Universal behavior was ob- 
served also for the hydrogen atom in a strong magnetic 



field pfj[, conformally deformed [|Tl|] and ray-splitting bil- 
liards |12l, and in the acoustic spectra of vibrating quartz 
blocks |13| ] . In all cases the general features of the conjec- 
tured universal behavior had been reproduced reasonably 
well, but a number of significant discrepancies remained 
unexplained. 

This was our motivation to study different types of 
billiard level dynamics a bit more detailed. All results 
to be presented below have been obtained in microwave 
billiards fl4|| . Here it is sufficient to note that for flat 
resonators the electromagnetic spectrum is completely 
equivalent to the quantum mechanical spectrum of the 
corresponding billiard, as long as one does not surpass 
the frequency v max = c/2h, where h is the resonator 
height. In the experiments we choose h = 8 mm yielding 
a maximum frequency of 18.74 GHz. 




FIG. 1. Sketch of the billiard used for the local level 
dynamics (in scale). 

One of the systems studied was a quarter Sinai billiard 
with a width b = 200 mm, a radius r — 70 mm of the 
quarter circle, and a length a which was varied between 
480 and 500 mm in steps of 0.2 mm. About 120 eigenval- 
ues entered into the data analysis in the frequency range 
14.5 to 15.5 GHz. The second system was a rectangu- 
lar billiard with side lengths a = 340 mm, b = 240 mm, 
containing 20 randomly distributed circular disks with a 
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diameter of 5 mm (see Fig. [j]). By a spatially resolved 
measurement |l5|] we found that all eigenfunctions tp in 
the studied frequency range were delocalized and \i/j\ 2 
was Porter- Thomas distributed. The position of one of 
the disks was varied in one direction in steps of 1 mm. 
Whereas the first type of level dynamics may be consid- 
ered as global, since a shift of the billiard length of the 
order of 1 wavelength will change the wavefunction pat- 
tern everywhere in the billiard, the shift of the disk gives 
rise to a local modification only. 




FIG. 2. Velocity distribution in a quarter Sinai bil- 
liard with one length as the level dynamics parameter. 
The solid lines correspond to a Gaussian distribution and 
a distribution described by a modified Bessel function 
(see Eq. (^)), respectively. The inset shows the distribu- 
tion in a logarithmic scale. 
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FIG. 3. Velocity autocorrelation function in a quarter 
Sinai billiard, where the level dynamics parameter was 
scaled according to Eq. (|l|). No less than 2000 velocity 
pairs entered every bin of the histogram. The solid line 
corresponds to the universal autocorrelation function of 
Simons and Althsuler. 

We start with a discussion of the global level dynamics. 
Figure § shows the velocity distribution for the quarter 
Sinai billiard with length a as the level dynamics param- 
eter. The distribution is well described by a Gaussian 
in accordance with the expected universal behavior (this 
result has been presented already in Q). Figure || shows 
the corresponding velocity correlator. To obtain the re- 
sult, each eigenvalue was studied over a range of four 



to five avoided crossings, and the scaling was performed 
by calculating the mean squared velocity for each eigen- 
value independently. Subsequently the results of about 
120 eigenvalues were superimposed. The solid line corre- 
sponds to Simons' and Altshuler's universal function Jl6| . 
The overall agreement between experiment and theory is 
good, but for x > 2.5 (not shown) the correlation func- 
tion does not approach zero but stays at negative values. 
This is an artifact resulting from an insufficient num- 
ber of data points making the calculation of the average 
( de n {X+x) g^m\ unreliable for large x values. Most 

correlation functions found in the literature end at x val- 
ues of at most 1.5, probably just for this reason. 
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FIG. 4. Velocity distributions in a rectangular bil- 
liard with randomly distributed disks with the position 
of one disk with diameter D as the level dynamics pa- 
rameter. The ranges of 5 = kD are 0.35 < S < 0.65 (a), 
1.4 < S < 2.6 (b), 5.1< S < 5.9 (c). 
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FIG. 5. Velocity autocorrelation functions in a rect- 
angular billiard with randomly distributed disks, where 
the level dynamics parameter was scaled according to 
Eq. (jl|). The ranges of S are the same as in Fig. ^j. 

Let us now turn to the discussion of the local level 
dynamics, where the position of one disk was varied. 
Whether a level dynamics must be considered as global 
or local, depends on the parameter S = kD, where D is 
the diameter of the disk, and k the wavenumber. It is 
well known that in the limit of small S values the spec- 
tral properties of billiards containing hard spheres devi- 
ate significantly from random matrix behavior \ \n\ . Fig- 
ure H shows the velocity distributions for three different S 
ranges. In Fig. |J(a) a disk with D = 5 mm was used, and 
the eigenvalues were taken in the frequency range 3.4 to 
6 GHz. In Figs. ||(b) and (c) the diameter of the movable 
disk was D = 20 mm with eigenvalues in the frequency 
ranges 3.4 to 6 GHz and 12.5 to 14.5 GHz, respectively. 
None of the found velocity distributions is Gaussian. One 
observes instead a distribution with a pronounced peak 



at v = 0, decreasing only exponentially for large values 
of \v\. With increasing 8 values the distributions turn 
gradually into a Gaussian. We completed the series by 
a level dynamics measurement for a half Sinai billiard, 
where the position of the half circle was varied. Here the 
obtained velocity distribution (not shown), correspond- 
ing to 5 values between 30 and 37, was already close to 
a Gaussian distribution. 

Figure |B| shows the corresponding velocity autocorre- 
lation functions. The scaling technique applied was the 
same as above. There is no longer any similarity between 
the experimental curves and the universal function. Only 
for the largest S value displayed, the experimental curve 
seems to approach the Simons- Altshuler correlation func- 
tion again. 

The results can be understood, if the movable disk is 
interpreted as a perturber probing the field in the res- 
onator (the perturbing bead method has been used many 
years ago to map the field distributions in microwave cav- 
ities and has recently been applied to the study of 
wavefunctions in chaotic billiards as well [fl9|~[?l"l| ) . In 
two-dimensional billiards the insertion of a metallic per- 
turber leads to a negative frequency shift proportional to 
E 2 , where E is the electric field strength in the resonator 
in the absence of the perturber. This holds as long as 
the dimensions of the perturber are small compared to 
the wavelength, i.e., in the limit 5 — * 0. Applied to the 
present problem this means that the eigenvalue velocity 
is given by dE n /dr — aV|-0| 2 where V is the gradient 
in the direction of the displacement, and a is a constant 
depending on the geometry of the perturber. It follows 
for the velocity distribution function 

P(v) = (5(v ~ 2cm/>VV0) • (3) 

Under the assumption that the wavefunctions can be 
described by a random superposition of plane waves [ p2[ , 
ib and are uncorrelated, and Gaussian distributed 
©, 

w) = V 2^~~ p ^ = H^-^. 

(4) 

The influence of the boundary is negligible here, since 
the linear dimensions of the billiard exceed the typical 
wavelength by factors of 5 to 10. Using Eq. (|]) the aver- 
age (||) is easily calculated and yields 

P(v) = ?-K o (0\v\), (5) 

where Ko(x) is a modified Bessel function, and f3 = 
A/2ak. The solid lines plotted in addition to the Gaus- 
sian curves in Figs. || and || have been calculated from 
Eq. (|J). In the limit of small 6 values distribution (|5|) 
describes the experimental distributions perfectly. 

The influence of local perturbations on the energy lev- 
els has been studied by Aleiner and Matveev p4] who 
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derived an explicit expression for the joint distribution 
function of initial and final energy levels. In their model 
the velocities are Porter-Thomas distributed pjj, if the 
coupling strength is taken as the level dynamics param- 
eter. The same distribution would have been expected 
in our case, if the coupling strength a would have been 
varied instead of the position (which, however, would be 
technically difficult to realize). 

For the quadratical average of the eigenvalue velocities 
we obtain using Eq. (S) 
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Entering with this expression into Eq. ([[]), we get for 
the rescaled parameter 
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where we have used that in billiards the mean level spac- 
ing is given by A = 4ir/A. Equation (0) shows that for 
the local level dynamics x is not an universal parameter, 
since it depends via a on the geometry of the movable 
disk. We shall therefore use the rescaled parameter 



kr 
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instead in the following. From the approach of random 
superposition of plane waves |2^| the velocity autocor- 
relation function can be easily calculated, too. Using 
standard techniques as they are described, e.g., in Ref. 
HI, we get 

c{x) = - [J 2 (5)]" = J 2 (x) - 2jf(x) - J (x)J 2 (x). (9) 




FIG. 6. Same as Fig. [](a), but with the level dynam- 
ics parameter scaled according to Eq. (^). 

Figure || shows again the velocity autocorrelation of 
Fig. @(a) for the local level dynamics, but now as a func- 
tion of x. The solid line corresponds to the theoretial 
expectation (Q). The experimental curve follows closely 
the predicted oscillations. With increasing 5 the oscilla- 
tions are more and more damped, but the wavelength is 
still in accordance with the theory (not shown). 



This paper has shown that two different regimes of 
level dynamics have to be discriminated. In the local 
regime velocity distributions and autocorrelation func- 
tions are quantitatively described by the approach of ran- 
dom superposition of plane wave, if the scaling (^|) is ap- 
plied. In the global regime, on the other hand, Simons' 
and Altshuler's universal functions describe the experi- 
mental results well, and the scaling (Q) is the appropri- 
ate one. The parameter <5 = kD governs the transition 
between the two regimes. 
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